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ABSTRACT 

A Monte Carlo code (artis) for modelling time-dependent three-dimensional spectral 
synthesis in chemically inhomogeneous models of Type la supernova ejecta is pre- 
sented. Following the propagation of 7-ray photons, emitted by the radioactive decay 
of the nucleosynthesis products, energy is deposited in the supernova ejecta and the 
radiative transfer problem is solved self-consistently, enforcing the constraint of energy 
conservation in the co-moving frame. Assuming a photoionisation dominated plasma, 
the equations of ionisation equilibrium are solved together with the thermal balance 
equation adopting an approximate treatment of excitation. Since we implement a fully 
general treatment of line formation, there are no free parameters to adjust. Thus a 
direct comparison between synthetic spectra and light curves, calculated from hydro- 
dynamic explosion models, and observations is feasible. The code is applied to the well 
known W7 explosion model and the results tested against other studies. Finally the 
effect of asymmetric ejecta on broad band light curves and spectra is illustrated using 
an elliptical toy model. 

Key words: radiative transfer - methods: numerical - supernovae: general. 



1 INTRODUCTION 

Type la supernovae (SNe la) are commonly believed to orig- 
inate from thermonu clear explosions of deg enerate material 
in white dwarf stars (|Hovle fc Fowleilll960l ). thus bemg im- 
portant for binary evolution and the chemical evolution of 
the Universe. By their dynamical interaction with the in- 
terstellar medium they may also play an important role in 
the star formation history of galaxies. However, even after 
many years of intensive research it is still not completely un- 
derstood how these explosio ns take place and which a re the 



ity and decay time-scale. iPhiUipsI (|l993l ) derived empirical 
relations to calibrate the peak luminosity of SNe la by dis- 
tance independent light curve properties which makes them 
standardiseable candles. Because of their high luminosities 
(typically ~ 10*^ ergs~^) they can thus be used to measure 
the cosmic expansion history up to high redshifts. Various 
forms of this approac h have been used by di fferent groups 
(e.g. IRiess et 



approac n nave Deen usea py gi rterent groups 
al.l ll998l : IPerlmutter et"ai1ll999l '). They found 
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progenitor systems [see e.g. iHillebrandt fc Niemeve^ (|2000|') 
for a review on the explosion mechanism and iLivid 
on progenitors]. 

In the standard single degenerate Chandrasekhar mass 
model, a C-l-O white dwarf accretes mass from a non- 
degenerate companion star (either a main-sequence or red 
giant star) until it nears the Chandrasekhar mass. Then the 
pressure of the degenerate electrons can no longer balance 
gravity and a thermonuclear runaway, disrupting the white 
dwarf, can occur. An attractive feature of this picture is 
that exploding at a well specified mass provides a natural 
explanation of the observed homogeneity of the SNe la class. 

However, SNe la are not all the same. Observations 
of nearby SNe la show a broad scatter in peak luminos- 
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that SNe la at high redshift appear fainter than expected 
in a matter dominated Universe, indicating an accelerated 
expansion. This was confirmed by following studies (e.g. 
iRiess et al.l 12004 lAstier et al1l2006h and led to the ACDM 
cosmologies in which a positive cosmological constant or 
"dark energy" with negative pressure is used to model the 
accelerated expansion of the Universe. 

One way to constrain dark energy models is to study 
the evolution of the dark energy equation-of-state param- 
eter w by measuring the c osmic expansion history , a goal 
of ongoing [e g. ESSENCE (|Wood-Vasev et al.ll2007l ). SNLS 
(|Astier et al.ll2006 ')] and future projects. For that purpose, 
the distances to high-redshift SNe la must be determined to 
high accuracy, which requires, amongst other factors (e.g. 
the extinction to the supernova and K-corrections) , a re- 
liable calibration of the peak luminosity. To improve the 
calibration techniques, which are based on purely empirical 
relationships, a thorough theoretical understanding of SNe 
la and improved observations are needed. 
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Much progress (both in numerical methods and compu- 
tational power) has been made in the hydrodynamical explo- 
sion modelling of SNe la since the first one-dimension a l (ID ) 
calculations by e.g. iNomoto. Thielemann fc Yokoil (1 19841 ) 
became available. Tod ay, fully three-dimensional (3D) 
explo si on models (e. g . iReinecke. Hi llebrand t fc Niemeveij 
20021: iGamezo et al.l l2003l : iRopke fc Hillcb randtl l2005l : 



Ropke fc Niemeve r^2007) are the state-of-the-art and have 



shown that 3D effects are essential to properly simulate 
the instabilities and turbulence effects which drive the ther- 
monuclear combustion. Furthermore they show that ejecta 
asymmetries can arise, either by hydrodynamical insta- 
biliti es during the burning phase or an asymmetric igni- 
tion l|Hoflich fc Stein|[2002l : lKuhlen. Wooslev fc Glatzmaie^ 
l2006l ). 

However these simulations - which give velocities, den- 
sities and composition of the explosion ejecta - are not di- 
rectly comparable to observations of real SNe la. For that 
purpose synthetic spectra and light curves must be obtained 
by radiative transfer calculations. This requires a solution of 
the multi-line transfer problem in expanding media where 
the opacity is dominated by the wealth of lines associated 
with the iron group elements which were synthesised in the 
thermonuclear explosion. Many ID studies have addressed 
this problem in the past either as s umin g pure resonance 
scatt ering (e.g. [Branch et al.l Il982l. 19831: Mazzali fc Lucvl 
1 19931 ) or p u re abso rption (e.g. iJefferv et al.l Il992l ') in the 
lines. ILucvI (Il999bl ) introduced an approximate treatment 
of line fluorescence. But that work is still in ID as are the 
studies done with t he general-purpose radiation transport 
code PHOEN I X (e.g. Lentz et al. 20011). Recently foUowing 
ILucvI (|2005l ). iKasen. Thomas fc Nugentl (|2006l ) described a 
3D time-dependent radiative transfer code which is capable 
of treating line fluorescence in an appro ximate way (similar 
to lLucvl[l999bl : |Pinto fc EastmanlbOOd ). 

Such studies have shown that to address the complex- 
ity of the hydrodynamic explosion models a 3D treatment of 
radiative transfer which simulates the 7-deposition and spec- 
trum formation in detail is needed. In particular, a careful 
treatment of the ionisation balance and a proper simulation 
of the redistribution of flux by li ne fluorescenc e [crucial for 
the near- infrared light curves, see iKasenI l|2006l )] is essential. 
Here we present a new Monte Carlo code (artis, Applied 
Radiative Transfer In Supernova e) which (based on the ap- 
proach of 'Lucv"200^, '2003', 2005') solves the time-dependent 
3D radiative transfer problem in chemically inhomogeneous 
models of supernova ejecta from first principles using a gen- 
eralised treatment of line formation and prioritising a de- 
tailed treatment of ionisation. The radiative transfer calcu- 
lation is parameter-free, depending only on the input model 
and atomic data, giving a maximum of predictive power for 
a given hydro model. Details of our code are given in Section 

m 

In Section[3]we use this code to calculate model spectra 
and light curves for the well-studied ID deflagration model 
W7 (|Nomoto et al.lll984l ) and compare them to observations 
and earlier synthetic results in order to test our code. We 
also investigate the influence of completeness of atomic data. 
Finally, in Section [J] line-of-sight dependent spectra and 
light curves for an ellipsoidal toy model are calculated to 
demonstrate the multi-dimensional capabilities of our code 



and study the effect of large scale asymmetries of the ejecta, 
before we draw conclusions in Section (5] 



2 METHOD 

Monte Carlo methods have been widely applied to model 
astrophysical radiative tr ansfer in ex panding media (e.g. 
lAbbott fc Lucy. ,19851 and ILucv fc Abbott 1993 for winds 



of massive stars and lMazzali fc Lucvl 1 19931 for the expand- 
ing envelopes of supernovae). This is mainly due to the 
ease with which they address the multi-line transfer prob- 
lem - a single photon can interact with many spectral lines 
due to the progressive redshift it experiences in the co- 
moving frame while it travels through the expanding en- 
velope. Monte Carlo methods are al so readily extende d to 
multi-dimensional problems (see e.g. iKasen et al. I I2OO6I and 
iMaeda. Mazzali fc Nomo to 2006) and easy to parallelise, al- 
lowing the exploitation of massively parallel super comput- 
ers. 

As the typical expansion velocities in supernova en- 
velopes (< 30 000 km s~^) imply large velocity gradients, the 
physical properties of the envelope will not change signifl- 
cantly within the resonance region of a spectral line (Doppler 
velocity vo ^ 5kms~^) and we adopt the Sobolev approx- 
imation and consider that line absorption of photons takes 
place only at a particular point of resonance. This simplifies 
the problem significantly and is a fundamental assumption 
of our study. For a detailed descrip tio n of the Sobolev ap- 

roxim ation see e.g. iMihalaj l| 19781 ) or lLamers fc Cassinellil 

19991 ). 

The second fundamental assumption we make is that 
of homologous expansion, thereby decoupling the radiative 
transfer from the hydrodynamic evolution of the ejecta. This 
applies if the ejecta are in free expansion, i.e. their kinetic 
energy density dominates the gravitational and internal en- 
ergy densities. For SNe l a this is achi eved less than a minute 
after the explosion (e.g. lRopkell2005l ). Thus it should be an 
excellent approximation for our radiative transfer calcula- 
tions which are usually started at ~ 2 days. In homologous 
expansion the velocity v of the ejecta at a particular posi- 
tion r is always proportional to the position, v — r/t with t 
being the time since explosion. 



2.1 Monte Carlo radiative transfer 

We extended the 3D Mo nte Carlo radiative transfer code 
introduced by ISimI (|2007l ) to a non-grey opacity treatment 
follow i ng th e schem e outlined in a series of papers by ILucvI 
(|2002l . |2003| . I2OO5I ). This scheme is based on the artificial 
division of the radiation field into indivisible energy packets 
as Monte Carlo quanta, rather than Nature's quantisation of 
radiation. This has two major advantages: first it keeps the 
code simple as packet histories can be followed one-by-one, 
avoiding the need to follow multiple packets in e.g. recom- 
bination cascades. Second it ensures rapid convergence to 
an accurate temperature stratification by imposing th e con- 
straint of energy conservation in the co-moving frame (|Lucy| 
Il999al ). 

In the following we outline the operation of the code 
with reference to the flow chart shown in Figure [1] and dis- 
cuss the relevant physical processes, giving detail only where 



Time- dependent 3D spectrum synthesis for type la supernovae 3 



we deviate from the formulations used bv lLucvl l|2002l , l2003l . 
I2OO5I ). The basic transfer calculations are carried out in the 
rest frame (rf, quantities always denoted unprimed) so that 
an energy packet's trajectory at time t is described by its po- 
sition r{t) and direction fJ.{t). For interactions with matter 
we transform to the local co-moving frame (cmf, quantities 
denoted with prime). 

2.1.1 Setting up the computational domain 

To avoid any symmetry assumptions and keep the code sim- 
ple, we map the explosion ejecta to a 3D Cartesian grid with 
origin at the centre-of-mass of the supernova and contain- 
ing 7'' cubic cells. The grid expands with time to follow 
the evolution of the ejecta. Physical variables are assumed 
not to vary spatially within the grid cells. As input we take 
densities Pi{t), velocities Vi{t) and composition Xi^k{t) from 
explosion models specified for the phase of homologous ex- 
pansion and map these for a time to onto the grid, to and all 
time specifications in the following indicate the time since 
explosion, unless otherwise noted. 

We then follow the expansion of the ejecta for A'' time 
steps up to time tjv by expanding the individual grid cells 
continuously. The time steps (f„,t„+i) are spaced logarith- 
mically and thermodynamic quantities in a cell such as den- 
sities, temperatures and atom/ion populations are kept fixed 
during a time step n at the value they have for tn+o.s,. 

2.1.2 Energy deposition 

SNe la light curves are primarily powered by the 
radioactive decays of '^ '^Ni -^^'^C o and ssgg^sep 
llTruran. Arnett fc CameronI Il967l : IColgate fc McKed 
1 19691 ) giving rise to the emission of a spectrum of 7-photons 
assoc iated with their transitions (|Ambwani fc Sutherland! 
Il988t ). Out of the total energies, Ewv and Eco emitted per 
decay of ^'^Ni and '^''Co, we determine the total 7-ray energy 
emitted in the decay chain for t ^ oo 

Etot = (^Ni -I- -Bco)A/Ni/mNi. (1) 

Mni is the initial mass of ^^Ni synthesised in the explosion 
and TTiNi the mass of the ^^ Ni atom. 

Following ILucvI (|2005l ). this energy is quantised into 
Af — -Etot/eo identical energy packets of cmf energy eo 
which we call "pellets" in the following. These pellets are 
distributed on the grid according to the initial ^^Ni dis- 
tribution and follow the homologous expansion until they 
decay. Decay times are sampled randomly according to the 
'^Ni— >^^Co— >^^Fe decay chain with pellets assigned to rep- 
resent either the ^^Ni or ^®Co decay. 

Upon decay, a pellet transforms to a single 7-packet rep- 
resenting a bundle of monochromatic 7-radiation with cmf 
energy eo and a cmf photon energy Ei.^ which is randomly 
sampled from the 7-lines in the appropriate decay of the se- 
quence '^^Ni^^'^Co ^'^'^Fe. The 7-packet's direction (/i') in 
the cmf is sampled randomly assuming isotropic emission. 

2.1.3 Propagation of ^-packets 

The 7-packets are propagated through the ejecta in the rest 
frame of the grid until either (i) they leave the grid, (ii) the 



current time step finishes or (iii) they interact with matter. 
Computing the distances to all possible events along the 
pac ket tra j ector y, we select the event which is reached first 
[see ILucvI (|2005l ) for details]. In the first case the 7-packet 
is fiagged inactive and the calculation proceeds to the next 
active packet. In the second case we save the rf data string 
(r, t, fi, e, E^) for the following time step and continue with 
the next active packet. 

For the third case we consider Compton scattering by 
free and bound electrons, photoelectric absorption and pair 
production as possible physical processes, the last being only 
available for 7-packets with energy E > 2m e(? . For details 
see again ILucvI (|2005l ) and Equation 1 of ISim fc Mazzalil 
( 2008 ) for the adopted phot oabsorpt ion cross-section. Which 
of the processes happens is determined randomly according 
to their absorption and scattering coefficients. 

In the case of a photoelectric absorption, the 7-packet 
energy is deposited as thermal kinetic energy. In the frame- 
work of this code, this is described by a transformation of 
the 7-packet into a so-called fc-packet of equal cmf energy. 
The treatment of fc-packets is described in section 1.41 

As our energy packets are indivisible, the treatment of 
Compton scattering and pair production, where the pho- 
ton energy is distributed to two particles, is slightly more 
complex. Following again ILucvI (|2005l , for Compton scat- 
tering the 7-packet is either scattered and continues as a 
7-packet of the same cmf energy as the incident packet or it 
is transformed into a non-thermal e~ -packet. e~ -packets are 
assumed to thermalise and are instantaneously transformed 
into fc-packets. For pair production we either create e^- or 
e~-packets. Assuming in situ annihilation, for an incident 
7-packet of cmf photon energy E'.y a fraction of 2mc(? / E'.^ 
(represented by the e'''-packets) is released in form of 7-rays 
at 0.511 MeV when the positron annihilates. The remain- 
der (representing the kinetic energy of the electrons and 
positrons) goes directly to the thermal pool. 

2.1.4 Treatment of thermal kinetic energy 

Neglecting energy storage in the ejecta gas, thermal ki- 
netic energy converts instantaneously (i.e. without propa- 
gating) into ultraviolet-optical-infrared ( UVOIR) radiation. 
This happens either directly via continuum emission by free- 
free or free-bound processes or indirectly by coUisional exci- 
tations/ionisations of the gas and subsequent radiative de- 
excitations/recombinations. In our framework, this means 
either transforming a fc-packet into an r-packet - represent- 
ing a monochromatic energy packet of UVOIR radiation (for 
details see Section [2. 1.6p - or into a packet of atomic inter- 
nal energy, which we call i-packets. Details of the latter will 
be discussed in the next section. 

Which of the processes happens depends on the cooling 
rate associated with the particular process and is sampled 
randomly. For the free-bound cooling rate we consider only 
spontaneous recombination since stimulated recombination 
is treated as negative photoionisation and, as such, is taken 
into account in the creation of fc-packets. 

From sampling the cooling rates, we know in which free- 
bound continuum (specified by level i of ion j of element 
k) to emit or whether free-free emission was selected. In 
both cases the r-packet's direction in the cmf fi' is deter- 
mined according to isotropic emission and the cmf energy 
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Figure 1. Flow chart outlining the mode of operation of the code. For discussion see text. 



e' of the generating packet is conserved. After transforming 
these quantities into the rf together with the packet's posi- 
tion (r), we must determine the r-packet photon frequency 
before continuing its propagation through the grid. 

For free-free emission the frequency is determined from 



10 



(2) 



with z a random number in ]0, 1[ and j^if) tx exp"'"'/''-^'' 
(we neglect the frequency dependence of the velocity- 
averaged Gaunt factor). For free-bound emission we use 
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(3) 



with z a random number in ]0, 1] and where Vi^j^k is the 
edge frequency of the continuum. The emissivity jf^j^^iv) 
for spontaneous r ecomb i nation by that continuum is given 
by Equation 23 of lLucvl (120031 ). 

In principle we want to only release thermal energy 
when transforming a fc-packet to an r-packet. However in 
emitting a bound-free photon at frequency u, a fraction 
i^i,j,k/v of the packet energy is drawn from the ionisa- 
tion/excitation energy pool. The machinery described in the 
next section ensures that this is balanced by the amount of 
thermal energy which is released by recombination pro cesses 
in ou r treatment of atomic internal energy (see also ILucvI 
l2003l ). 



2.1.5 Treatment of atomic internal energy 

We treat atom i c internal e nergy using the macro-atom ap- 
proach of lLucvl (|2002l . l2003l ) which is based on an interpreta- 
tion of statistical equilibrium in terms of macroscopic energy 
flow rates. For Monte Carlo metho ds this is an improve- 
ment over e arlier wo rk in ID (e.g. iMazzali fc Lucvl 1 19931 : 
lLucvlll999bl ) or in 3D (|Kasen et aLll2006l '). We stress that the 
macro-atom approach allows for a fully general treatment of 
radiation/matter interactions in statistical equilibrium, in- 
cluding true absorption, scattering and fluorescence. 



Whenever a packet is converted to atomic intern al en- 
ergy (z-packet), we assign it to a macro-atom state (see lLucvl 
|2002| for relevant definitions). The ejecta gas can be ex- 
cited/ionised either coUisionally, represented by a fe-packet 
to z-packet transition (see last section) or radiatively by 
UVOIR radiation, represented by an r-packet to i-packet 
transition. Non-thermal excitations by fast electrons and 7- 
rays are not currently taken into account but the method 
could be extended to include them in due course. 

For fc-packet to i-packet transitions, the macro-atom 
state is determined by sampling the cooling rates for coUi- 
sional excitation/ionisation, which are c alculated in th e van 
Regemorter and Seaton approximations l|Mihalasll 19781 ). For 
r-packet to i-packet transitions, the absorption process iden- 
tifies the appropriate macro-atom state directly. Details of 
the treatment of r-packets are given in Section [2. 1.61 

We do not take into account energy storage in the ejecta 
gas but assume i-packets convert instantaneously back into 
k- or r-packets. For that purpose we calculate rates for all 
transitions which connect an activated state to other macro- 
atom states or which allow the macro-atom to deactivate, i.e. 
transform the internal atomic energy back to either thermal 
kinetic (fc-packet) or radiative (r-packet) energy. From all 
these transitions we randomly select one by samplin g the 
probabilities which are computed following ILucvI (|2003l ). If it 
is a macro-atom internal jump we reset the i-packet macro- 
atom state and continue this procedure until we select a 
deactivating process. 

During deactivations we conserve the cmf energy of the 
packet to enforce radiative equilibrium. In the case of a col- 
lisional deactivation the i-packet is just transformed into a 
fc-packet. In case of radiative deexcitation/recombination we 
transform the i-packet into an r-packet and emit it isotrop- 
ically in the cmf. Finally we need to assign a photon fre- 
quency to the r-packet. For a radiative deexcitation this is 
given by the frequency of the associated line transition in 
the cmf. In the case of a radiative recombination, we again 
sample the emissivity for the free-bound emission (Equation 
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[3|. Note that we thereby emit both ionisation and thermal 
energy, as we do it in the case of free-bound cooling. 

2.1.6 Propagation of UVOIR radiation 

r-packets representing monochromatic energy packets of 
UVOIR radiation are propagated through the ejecta from 
their first emission by any of the above processes until (i) 
they leave the grid, (ii) the current time step finishes or 
(iii) they interact with matter. In the zeroth time step we 
also have radiation which comes from pellets which decayed 
before to. We assume that the ejecta are optically thick at 
these early times such that the radiation is trapped and 
local thermodynamic equilibrium (LTE) is a good approxi- 
mation. Therefore we emit these packets isotropically at the 
position the pellets would have at to and assign them fre- 
quencies according to a black-body distribution at the local 
kinetic temperature Tc. 

We follow the packets along their trajectory until they 
are stopped by the first event which occurs. For (i) and 

(ii) this is a simple geometrical problem. For case (iii) it 
requires calculating the distance to a randomly sampled 
optical depth Ti = — In z with z €]0, 1], as described by 
iMazzaU Lucvl (|l993l l. First the trajectory point at which 
the photon comes into resonance with the next spectral line 
of the line list is calculated in the Sobolev approximation. If 
the continuum optical depth accumulated up to that point 
is bigger than the random optical depth, the travelled dis- 
tance ds is calculated from the continuum opacity (Kcont) 
from Kcontds = Tr and a continuum absorption occurs. If 
the sum of continuum and line optical depth is bigger than 
the random optical depth, ds is given by the distance to the 
resonance point and a line event occurs. Otherwise we cal- 
culate the distance to the next spectral line with which the 
photon comes into resonance and repeat the procedure. 

According to the chosen event, the packet properties are 
now updated. For case (i), the r-packet is fiagged inactive 
and we proceed to the next active packet. For case (ii), we 
save the actual rf data string (r, t, fi, e, u) for the following 
time step and continue with the next active packet. For case 

(iii) , we change the packet type according to the process 
which terminated the packet flight. For a line absorption, 
we convert the r-packet into an i-packet and set it to the 
corresponding macro-atom target state and the macro-atom 
machinery (Section 12. 1.5|l determines its fate. 

For a continuum event, we first decide whether it was 
Thomson scattering, free-free absorption or bound-free ab- 
sorption by randomly sampling their absorption coefficients. 
Thomson scattering, which we assume to be coherent and 
isotropic, changes neither the packet type, nor the cmf fre- 
quency of the packet. Thus, Thomson scattering amounts to 
isotropic reemission of the r-packet in the local cmf. Free- free 
absorption transforms the absorbed radiative energy com- 
pletely into thermal kinetic energy, i.e. converts the r-packet 
into a fc-packet. 

Bound-free absorption is more complicated since it con- 
tributes both to atomic internal energy and thermal kinetic 
energy. After randomly selecting one specific continuum ac- 
cording to its absorption coefficient, we follow Lucy in forc- 
ing a conversion to either a single k- or i-packet. For an 
r-packet absorbed with frequency v hy & continuum with 
edge frequency i^i,j,k, a fraction Vij^k/i^ of its energy is con- 



verted to atomic internal energy and the rest to thermal 
kinetic energy. To account for this, we convert the packet to 
an j-packet with probability Ui_j_k/v, setting it to the corre- 
sponding macro-atom state, and to a fc-packet otherwise. 

To speed up the calculations, we introduced an op- 
tional initial grey approximation in the code. Controlled by 
two input parameters TTh.min and Ngrey, the propagation 
of UVOIR radiation in cells which have a Thomson opti- 
cal depth TTh > TTh.min is performed including Thomson 
scattering as the only source of opacity. If tth falls below 
TTh.min, aud for all time steps later than Ng^oy, we switch 
back to the non-grey treatment. This speeds up the initial 
phase of a calculation by a factor of ~ 10, since much of the 
ejecta are still rather dense and photons cannot propagate. 
With a good selection of rTh,min and Ngi^y the differences 
in spectra and light curves from a fully non-grey calculation 
are negligible. 



2.1.7 Extraction of spectra and light curves 

When the simulation has finished, we extract the spectral 
evolution by binning the escaping r-packets in frequency, 
time and direction of escape, accounting for light travel-time 
effects. Colour light curves are extracted from the spectral 
evolution by i ntegrating the spectra over the ap propriate 
filter functions l|Bessell fc Brett|[l988l : lBesselllll990l l. UVOIR 
bolometric light curves are extracted by binning the escaping 
r-packets by time and angle. Similarly 7-ray spectra and 

light curves can be obtained from the 7-packets. 

Us ing the formal integral method developed by ILucvI 
(|l999bl ) to extract spectra and light curves would substan- 
tially reduce the Monte Carlo noise. However, the need to 
store line and continuum source functions in all grid cells for 
all lines and continua included in the simulation (typically 
10^ . . 10^) causes the simulation memory requirements to 
be so high that this method becomes computationally pro- 
hibitive, even with current state-of-the-art super computers 
for 3D grids. Therefore all results presented here have been 
obtained using the direct binning approach. 



2.2 Plasma conditions 

To calculate opacities we need to know atomic level popula- 
tions which in turn depend on the radiation field. In princi- 
ple, these could be extracted exactly out of our simulation, 
but it is computationally restrictive to record the complete 
set of level-by-level radiative rates for every grid cell which 
would be required to do this. Therefore we use an approx- 
imate model for the radiation field and for excitation and 
ionisation conditions. Two different descriptions have been 
implemented in our code. The first parameterises the ionisa- 
tion balance in a simple one-temperature description using 
the Saha formula. This is less demanding of computational 
resources but also less accurate than our second which treats 
ionisation in detail and simultaneously solves the thermal 
balance equation. 
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2.2.1 Radiation field models 

Following ILucvI l|2003l ) , the zeroth moment of the radiation 
field in a grid cell can be reconstructed from the Monte Carlo 
packets using the estimator 

1 



AT,. 



Ecmf J 
e„ As 



(4) 



where At is the length of the time step, V the grid cell 
volume, ds the trajectory length within the grid cell and e""' 
the packet energy in the cmf. The summation runs over all 
trajectory segments ds in V for which the packet frequency 
v & [u,v + Au\. 

As discussed above, we do not use this general descrip- 
tion directly but parameterise the radiation field. For the 
simple one-temperature description of ionisation, the radia- 
tion field is parameterised as a black body. 



J. = (Tj) 



(5) 



at a temperature Tj corresponding to the local energy den- 
sity of the radiation field {B^{T) is the Planck function). 
The parameter Tj can be extracted by equating the fre- 
quency integrated estimator (J) = ,Ju Au (where ,Ju is 
given by Equation |4]) with the Stefan-Boltzmann law 



T, = (tt (J) /af'^ . 



(6) 



In our simple ionisation treatment, we also equate the kinetic 
temperature To = Tj. 

For the detailed ionisation treatment we parameterise 
the radiation field using a nebular approximation 



J. = WB, (Tr) 



(7) 



with the radiation temperature Tr and the dilution fac- 
tor W as parameters. Tr is chosen such that the mean- 
intensity weighted mean frequency of the radiation field 
(y) = Jg°° vJv A v I fp°° Jij Av matches th at of a blackbody 
at Tr. Following iMazzali fc Lucvl (119931 ) Tr and W can be 
extracted from (J) and {y) by 



W = 



and 



Tr^ 



xkB 



(8) 



with X = 360-(^ (5) /tt* and C, (x) being Riemann's (^-function. 
Tj remains defined as in the simple ionisation treatment 
(implying Tj = W'^^^Tk from Equation El . 



2.2.2 Excitation and ionisation 

Independent of the ionisation description used, we apply the 
Boltzmann formula evaluated at Ti 



^i,j,k 



(9) 



— ^-^^ exp — — 

no,j,k 90,j,k \ kbTi 

to calculate the population Uij^k of a level i of ion j of 
element k relative to the ion ground state population no,j,k. 
gij^k and eij-,fc are statistical weights and energies of the 
corresponding levels. 

Assuming ionisation equilibrium in a radiation domi- 
nated environment we require 



^ Currently wc consider only photoionisation and radiative re- 
combination. Temperatures in tlie recombination zones of iron 



(10) 



where Afj.k is the number of levels associated with ion j of 
element k, and rio the density of free electrons. "/i,j,k denotes 
the photoionisation rate coefficient from a bound level i,j, k 
and al'^j ^ is the rate coefficient for spontaneous recombina- 
tion fi;om__tlie_(j(-|- l)-th ion to the bound state i,j,k (see 
e.g. lMihalaj|l978l ). In our simple one-temperature ionisation 
treatment this reduces to the Saha equation 



J,k 



Nj+i^kne 
where 



Uj,k C 



exp 



A:bTc 



2 \27rmefcB 



3/2 



(11) 



(12) 



and Uj^k is the partition function of ion j of element k. 

For the detailed ionisation treatment we neglect stim- 
ulated recombination, which is small compared to sponta- 
neous recombination, and use the partition function together 
with 



- 3,k 



j: 



^i,j,k^i,j,k 
Tl'0,3,k 



and 



4sp 



E 



(13) 



to obtain the following ionisation formula from Equation 1101 



N,,k 



/Isp 



Uj,k 



Nj + l^kTlc Tj^k 90,j,k 



(14) 



Since a^^^ j, does not depend directly on the radiation field, 
we calculate Aj^k by summing all the a^J^ j, of the ion j, k 
evaluated at the local kinetic temperature T^. 

In contrast, the photoionisation rate coefficients 7ij,fc 
depend directly on the radiation field and should be deter- 
mined from volume-based Monte Carlo estimators for most 
accurate results. However, we cannot afford to store 7i,j,fe 
for each bound-free continuum in every grid cell. Instead we 
derive 7*j j(. from an integration of the radiation field model 
(Equation[7)l using the local values of Tr and W . Estimated 
in this way 7*j^fc is expected to be reliable for bound-free 
absorptions which lie around the peak of the radiation field. 
However, the radiation field model is inadequate at very blue 
frequencies where the spectrum is systematically affected by 
bound-free edges. Thus to track the effect of absorption con- 
tinua on the radiation field, we record estimators 7o,j,fc for 
the photoionisation rate coefficient of the ground level con- 
tinua and use them to derive renormalisation coefficients 
Co,i,fc = 7o,j,fc/7o,i,fc- To obtain the 7ij-,fc in Fj.fc of Equation 
I13l we renormalise all the integrated photoionisation rate co- 
efficients via 



> * 

li,j,k — ',0,j',fc'7i,i,fe> 



(15) 



group elements are sufficiently low that dielectronic recom bi- 
nation (calculated according to IShull fc van Steenberd [iQSd) is 
small compared to radiative recombination. However we note that 
in high temperature environments it may become important, es- 
pecially for intermediate mass elements. For more detailed studies 
our treatment could readily be extended to incorporate additional 
ionisation/recombination processes (including dielectronic recom- 
bination and non-thermal ionisation). 



Time- dependent 3D spectrum synthesis for type la supernovae 7 



where we take the (^o,i',fc'-value from the ground level con- 
tinuum which lies closest in edge frequency to the considered 
continuum (i.e., the one for level i,j, k). For continua redder 
than the reddest ground level continuum, ^ is set to 1. 

For the zeroth time step, and also for cells treated in 
our initial grey approximation (see Section I2.1.6p . we use 
the simple ionisation treatment. We also assume that at all 
times prior to our zeroth time step the ejecta are sufficiently 
optically thick, that all the radiation is trapped. Thus we 
calculate initial temperatures T^.o = Tr.o = T!j,o for each 
grid cell from the energy released by the decay of ^^Ni in 
the cell prior to the midpoint of the zeroth time step. Since 
the trapping is not really perfect, it takes a few (typically 
~ 5) time steps until the propagation of the energy packets 
washes out the effect of this boundary condition on temper- 
atures and level/ion populations. 

In practise this boundary condition primarily effects cal- 
culations using the detailed ionisation description. In this 
case, it is more important to have accurate level populations 
as errors can feed back to the thermal balance calculation. 
Thus we typically do the first ~ 10 time steps of a calcu- 
lation in the simple ionisation treatment to obtain reliable 
initial conditions for the detailed ionisation treatment at all 
later times. Since our simple ionisation treatment reproduces 
LTE conditions we expect it to be a good approximation at 
early times when optical depths are high. 



2.2.3 Thermal balance 

In the detailed ionisation description, we also solve for the 
local kinetic temperature Te in a grid cell by balancing the 
local heating TC and cooling C rates. The heating rates are 
calculated from the previous time step using volume-based 
Mont e Carlo estimators for the heat ing by 7-ray s (see ILucvI 
l2005h and free-free absorptions (see ILucvI|2003|) and event- 
based Monte Carlo estimators for the heating by collisional 
deexcitation / recombination. 

We cannot afford to store a bound-free heating estima- 
tor for each continuum in every grid cell. Thus bound-free 
heating rates 



'^i,3,k {if.j.k - li,j,k^ hVi 



(16) 



(17) 



where we again take the Co,i',fc' which lies closest in fre- 
quency to the considered continuum (i.e., the one for level 
i,j,k) to track the effect of absorption continua on the ra- 
diation field. Continua which are redder than the reddest 
ground-level continuum are not renormalised. 

We include cooling by free-free and free-bound emis- 
sion, collisional excitation/ionisation and adiabatic expan- 
sion. Since these cooling rates do not depend on the radia- 



tion field directly, they are calculated fro m the current set 
of level populations and Te as described in lLucvl l|2003l '). The 
adiabatic cooling rate is given by 



dV 
'IT 



(18) 



(with "yfj^k th e mod i fied r ate coefficient for photoionisation 
as defined by ILucvI [200^ are not obtained directly from 
Monte Carlo estimators but from an integration of the radi- 
ation field model using the local values of Tr and W yield- 
ing hYj\. As for the photoionisation rate coefficients (Sec- 
tion [2]2]2]), we then use estimators /ioj_fe for the bound- free 
heating coefficient of the ground level continua to derive 
renormalisation coefficients ^oj.fc ~ f^oj,k/h^^j*k- These are |200 ' 
used to calculate the bound-free heating coefficients 



with p the gas pressure and V the volume of the grid cell. 

Balancing the heating and cooling rates on a temper- 
ature interval [Tmin, Tmax] with Tmin — 3500 K and Tmax = 
1.4 • 10^ K in our simulations, leads to a new value of Tc. For 
cells in which the heating and cooling rates cannot be bal- 
anced on this interval. To is set to the upper/lower boundary 
value if the heating/cooling rates dominate. The same tem- 
perature limits apply to Tj and Tr. Tmax is chosen such that 
it will be only reached at the very earliest times in the inner 
iron-rich core where the opacity is so high that the radiation 
is trapped (actually for the simulations described in this pa- 
per Tmax was never reached). Tmin will be reached at late 
times but only in layers which are sufficiently optically thin 
that they have no strong infiuence on the spectrum forma- 
tion. 

Assuming that time steps are short enough and the ra- 
diation field does not depend strongly on the evolution of 
Tc we use the obtained value for Te and the same Monte 
Carlo estimators as before to re-solve for the populations and 
then for thermal balance again. Repeating this iteratively we 
reach a converged solution for To and the ionisation state. 
This iteration can be extended to an out er loop in wh i ch the 
Monte Carlo simulation is also repeated. iKasen et all (|2006l . 
Section 3.4) showed that, even from crude initial conditions, 
such an iteration converges very rapidly. By using short time 
steps and initial thermodynamic quantities from the previ- 
ous time step, we found that such an outer iteration of the 
Monte Carlo experiment does not affect the results strongly. 
Thus we do not iterate on the Monte Carlo simulation for 
each time step in the simulations described below (however 
the option to do so remains in the code). 



2.3 Model atoms 

The most important opacity source in supernova envelopes is 
the wealth of bound-bound line transitions associated with 
the iron group elements synthesised by the thermonuclear 
explosion. Thus an extensive line list which covers all the 
important species is needed. We restrict our atomic data to 
the lowest five ionisation stages of elements up to Zn, ne- 
glecting H, Li, Be and B. To study the sensitivity to the 
choice of atomic data , we use two da t asets which are based 
on the ~ 5-10^ lines of lKurucz fc Belli (|l995l l (CD23) and the 
more modern theoretical atomic data computed bv lKuru"c3 
20060 which are far more comprehensive but still not com- 
plete. For example this dataset contains ~ 3.5 • 10^ lines as- 
sociated with the important second and third ions of Fe, Co 
and Ni compared to only ~ 1.2 ■ lO'^ in CD23. 

Since it is computationally too demanding to use all 
of the lines in a dataset, when constructing model atoms 
we apply cuts in log (gf) to reduce the number of lines. 
From the comprehensive dataset we take only the lines of 
Fell, Felll, Con, Com, Nill and Nilll, the most important 



Available at http : //kurucz . harvard . edu/atoms . html . 
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Table 1. Model ato ms used i n the calculat ions. Source A is 
iKurucz fc Beiil ^99^), Source B iKuruc j 1120061 ). 



Model name 


Source 


log(9/) 


Levels 


Lines (total) 


cd23^f-2 


A 


-2 


18237 


1.4 ■ 10^ 


cd23^f-3 


A 


-3 


18815 


2.5 • 10^ 


cd23^f-5 


A 


-5 


19472 


4.1 ■ 10^ 


cd23.gf-20 


A 


-20 


21100 


4.7 ■ 10^ 


big_gf-3 


B (Feii-iii) 


-3 


41829 


3.5 ■ 10** 




B (Coii-iii) 


-3 


31816 






B (Niii-iii) 


-3 


24853 






A (other) 


-3 


16616 




big-gf-4 


B (Fell-Ill) 


-4 


47141 


8.2 ■ 10*^ 




B (Coii-iii) 


-4 


36442 






B (Niii-iii) 


-4 


28022 






A (other) 


-4 


17523 





species in the spectrum forming region, and continue to use 
the data of CD23 for all the other species and remaining ions 
of Fe, Co and Ni. In this way, we created a set of model atoms 
(summarised in Table [TJ which are used in the following 
calculations. 

Bound-free cross section s are obtained f rom the fits 
by IVerner fc YakovlevI (Il995l ) or IVerner et al] ([l996) where 
available. A hydrogenic approximation is adopted for excited 
configurations. 



3 APPLICATION TO THE W7 MODEL 

As a first test we calculate the spectral evolution for 
the one-dimensional deflagration model W 7 l|Nomoto et al.l 
1 19841 : [Thielemann. Nomoto fc Yokoil Il986l ) which has al- 
ready been investigated i n several other radi ative transfer 
studies (e.g. IJefferv et al.|[l992: Hoflich 1 995: ,Nugent et al.l 
19971: iLentz et al.ll200ll: ISalvo et aI]|200li:lBajon et al.ll200l 



Kasen et al.l 20061 ) and is found to be in good overall agree- 



ment with observations. 

Our calculations use 5 • 10^ energy packets to follow the 
spectral evolution over 100 time steps from 2 to 80 days after 
explosion (Alog(t) ~ 0.037) and were performed on a 50'^ 
grid using the cd23_gf-5 model atom (see Table [l]for details) 
with our detailed ionisation treatment. The first 10 time 
steps have been calculated in the simple ionisation treatment 
to obtain reliable initial values (see the discussion in Section 
I2.2.2|l . To save computational time we make use of our initial 
grey approximation setting TTh.min = 15 and A'grey = 40. 

Figure [2] shows the 2500 to 10 000 A spectra for 15, 22 
and 44 days after explosion, corresponding to -5, +2 and 
-|-24 days relative to B band maximum. The colour cod- 
ing shows the element associated with the last line emission 
of escaping photons, giving an indication of the elements 
responsible for spectral features (absorption features of in- 
dividually strong lines can be identified in the plot from 
their associated P Cyg emission). The contribution of pho- 
tons which escaped after a bound-free or free-free emission 
to the total flux is negligible (< 10""^) and not shown. The 
most striking individual features are the Ca ll near- infrared 
(NIR) triplet at ~ 8600 A and the characteristic Si ll line at 
6355 A which is clearly visible around maximum light. While 
the outer layers, which are dominated by intermediate mass 
elements, are optically thick, the spectra are dominated by 



atomic number 



7 6 7 8 10 12 14 16 18 20 22 24 26 28 30 

Day -5 




4000 6000 8000 

Wavelength (A) 

Figure 2. Spectra for the W7 explosion model at 5 days before 
and 2 and 24 days after maximum light in B band (from top 
to bottom) calculated using the detailed ionisation description 
and atomic dataset cd23_gf-5. Overplotted in blue are the ob- 
served spectra of SN 1994D for comparison (no redshift and ex- 
tinction corrections have been applied). The colour coding shows 
the element associated with the last line emission of escaping pho- 
tons. Atomic numbers in the colour legend are associated with the 
colour right to the number centre. 



these elements. Around maximum light the outer layers be- 
come optically thin and the spectra start to be dominated by 
Fe group elements which completely take over at later times. 
Only individual strong lines of the lower mass elements (e.g. 
the Call NIR triplet) persist. 

Overplotted w ith our synthetic spectra are observa- 
tions of SN 1994D (jPatat et al.ll 19961 ) for the corresponding 
epochs. Given that the W7 model has not been tuned to 
any particular supernova, the agreement of our model spec- 
tra with the observed ones is good. We reproduce the main 
spectral features, e.g. the Sill line and the Can NIR triplet, 
as well as the overall flux distribution. However there are 
obvious differences: at early times our model spectra have 
a strong excess in the UV fiux below 4000 A where Fen, 
Felll, Com and Coiv dominate the spectrum formation. At 
late times we obtain emission in Fe ll at ~ 6050 A where 
the data shows an absorption featu re. Model spectra calcu- 
lated for W7 byll asen et al] (|2006l ) show the same discrep- 
ancies compared to SN 1994D. This could be an indication 
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Figure 3. The left panels show the radial ionisation structure of Fe in the W7 explosion model at 31 days after the explosion obtained 
with the simple (top) and detailed ionisation treatment (bottom) and model atom cd23_gf-5. Black/red/green/blue/yellow lines represent 
Fe l/ll/lll/iv/v. The outer layers of the model, which consist only of unburncd material and extend up to 22 800kms~^, are not shown. 
The right panels show the time evolution of the Fe ionisation structure at a velocity of 9590 km s~^ (indicated by the grey vertical line 
in the left panels) from 2 to 70 days for the two different ionisation treatments. 



that some details of the explosion model would need to be 
changed to obtain better agreement with SN 1994D. 



3.1 Simple versus detailed ionisation treatment 

In the following we investigate how the two different ionisa- 
tion treatments (see Section I2.2.2p which are implemented 
in our code affect the spectral evolution of a given explosion 
model. Therefore we also followed the spectral evolution of 
the W7 model using the cd23_gf-5 atomic data from 2 to 80 
days after the explosion using the simple ionisation treat- 
ment. 

Since differences in the spectral evolution must orig- 
inate from the different ionisation treatment, we compare 
ionisation fractions from our different simulations in Fig- 
ure [3] We focus on Fe since the iron group elements, which 
have similar ionisation structure to each other, dominate the 
spectra for the relevant epochs. The left panels of Figure [3] 
show the ionisation fractions as a function of radial velocity 
at 31 days after the explosion for the simple and detailed 
ionisation treatment. The right panels show the same ion- 
isation fractions as a function of time at a radial velocity 
of 9590 km s~^, which is marked by the grey vertical line in 
the left panels. This is the velocity at which the iron group 
mass fraction of the ejecta drops to 0.5 and thus is around 
the outer velocity of the iron-rich inner core. 

From the left panels we see that, using the detailed ioni- 
sation treatment, the ejecta are more highly ionised at higher 
radial velocities. The right panels show that, with the de- 
tailed ionisation treatment, the ejecta also stay more highly 
ionised for a longer time at a given radius and never re- 



combine as fully as with the simple ionisation treatment at 
late times. In the simple ionisation treatment, the ionisation 
is determined only by the energy density of the radiation 
field calculated from the r-packets. In the detailed ionisation 
treatment, the frequency distribution of the calculated radi- 
ation field has an important effect. Tr > Tj and W < 1 (see 
Equations 6-8) with Tr and W departing more strongly 
from Tj and 1, respectively, with increasing time and de- 
creasing ejecta density and opacity. Thus, there tend to be 
more high-energy ionising photons than suggested by the ra- 
diation field model of the simple ionisation treatment. This 
gives rise to higher ionisation states, particularly for later 
times. This is illustrated in Figure U which shows the time 
evolution of the radiation temperature Tr and the tempera- 
ture Tj corresponding to the energy density of the radiation 
field at a radial velocity of 9590 km s~^ . While Tj drops simi- 
larly for the simple and detailed ionisation treatment, Tr for 
the detailed ionisation treatment stays significantly hotter 
from 20 days on. 

To shows a similar decline and stays somewhat below Tr 
up to about maximum light (note that at very early times, 
prior to ^ 10 days, this region is dominated by Fe v for which 
line transitions were not included - this likely explains the 
less complete coupling of Tc to the radiation field at these 
early epochs) . Afterwards 7-ray heating dominates the total 
heating rate and controls the kinetic temperature. Although 
the 7-ray heating decreases smoothly with time, Te shows 
distinct changes. These arise when recombination from one 
dominant ion to the next occurs and the contributing cooling 
processes change abruptly. In order to retain the balance 
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Figure 4. Temperature evolution at a radial velocity of 
9590 km i.e. at the outer edge of the iron-rich core. The green 
crosses show Tj in the simple ionisation treatment. The circles 
are the three temperatures used in the detailed ionisation treat- 
ment and show the kinetic temperature Te (black), the radiation 
temperature (blue) and Tj (red). At around 45 days, Tj hits 
the lower temperature boundary of 3500 K and is not allowed to 
drop further. 
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Figure 5. Radial temperature distribution at 31 days after the 
explosion. The outer layers of the model, which consist only 
of unburned material and extend up to 22 800 km s~^, are not 
shown. The green crosses show Tj in the simple ionisation treat- 
ment. The circles are the three temperatures used in the de- 
tailed ionisation treatment and show the kinetic temperature Tc 
in black, the radiation temperature T^i in blue and Tj in red. 
The black/violet/blue/green/red/orange vertical lines show the 
mean radii of last scattering in bolometric/CZ/B/V/i?// light. 
The sharp drop in Tc around 12 500 km s~^ is associated with a 
strong change in composition in the W7 model. The iron group 
abundance abruptly drops at this point and intermediate-mass 
elements dominate. 



between the heating and cooling rates then Tc must change 
rapidly. 

Figure [5] illustrates the radial temperature distribu- 
tion at 31 days after the explosion. Tj for the two differ- 
ent ionisation treatments shows no strong difference. How- 
ever Tr in the detailed ionisation treatment differs strongly 
from Tj over the entire ejecta, indicating that the radia- 
tion field is bluer and more dilute that would be implied 
by our radiation field model in the simple ionisation treat- 
ment. The kinetic temperature is mainly controlled by 7-ray 
heating which dominates over the other heating rates out- 
side ~ 2500 km s~^ and varies only gradually with velocity, 



declining outwards. This leads to Te > Tr. Around recom- 
bination fronts the contributing cooling processes change 
abruptly. This means that Tc has to change rapidly in such 
regions in order that the heating rate remains balanced by 
the total cooling rate and leads to the rapid decline in To be- 
tween 8000 and 10 000 km/s where the recombination from 
Felll to Fell happens (compare Figure [3)|. The vertical lines 
show the mean radii of last scattering for packets which had 
their last scattering at ~ 31 days. Because of line blocking 
the optical depths are biggest in the UV and blue and the dif- 
ferent bands probe deeper into the ejecta in a sequence from 
blue to red. An exception is the mean radius of last scatter- 
ing in the / band which lies at ^ 8500 km s~^, between the 
R and V bands. This is because the / band has significant 
contribution from the Ca ll NIR triplet which forms outside 
the iron-rich core. 

That the ejecta stay more ionised for longer in the de- 
tailed ionisation description has a direct influence on the 
spectra at late times. This is illustrated in Figure [6] which 
shows spectra for 15, 20, 40 and 60 days after the explo- 
sion for the simple and detailed ionisation treatments. While 
there are no differences up to about maximum fight in B 
band at 20 days, the spectra calculated using the detailed 
ionisation treatment are significantly bluer at late times. 
This is because they still have more Fell contribution rela- 
tive to Fe I than the spectra using the simple ionisation treat- 
ment (Fell has bluer lines). The simple ionisation treatment 
shifts more flux into the red and NIR leading to a strong Fe I 
emission feature at 60 days which is not seen in observations. 

The same effect influences the broad band light curves 
calculated using the simple and detailed ionisation treat- 
ment which are s hown in Figure [T] Th e band passes of 
iBesselll (|199G| ) and iBessefi fc Brett! lfl988l ) are shown in the 
top panel of Figure |S]to clarify the following discussion. Like 
the spectra, the light curves are rather similar before maxi- 
mum hght in B band at ~ 20 days, with slight differences in 
the NIR bands. After maximum light there are strong dif- 
ferences. The U , B and V fight curves calculated using the 
simple ionisation treatment fade much quicker than those 
using the detailed ionisation treatment. In contrast, R, I, J, 
H and K stay brighter with J showing a third maximum 
which is associated with the appearance of the strong Fel 
emission feature in the bottom right panel of Figure [§] at 
about 60 days. 

Of particular interest is the R band which shows a clear 
secondary maximum using the detailed ionisation treatment 
but only a slight platea u using the sim ple ionisation treat- 
ment. As pointed out by iKasenI l|2006l ). the secondary max- 
imum forms when the zone in which doubly ionised iron 
group elements recombine hits the inner iron-rich core. It is 
then that the redistribution of flux from the UV and blue 
part of the spectrum into the red and NIR by fluorescence 
is most effective. This argumentation is confirmed and ex- 
plains the differences in the R band light curve between the 
simple and the detailed ionisation treatment. In the right 
panels of Figure |3l which shows the ionisation fractions of 
Fe as a function of time at the outer edge of the iron-rich 
inner core, we see that for the simple ionisation treatment 
(top panel) the transition from Fe ill to Fe ll happens at ~ 33 
days while it occurs at ~ 36 days in the detailed ionisation 
treatment. These times correspond to the times of the sec- 
ondary maximum of the R band light curve in Figure [T] 
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Figure 6. Spectra for the W7 explosion model at 15, 20, 40 and 60 days after the explosion, respectively. The green and red lines 
were calculated with the cd23_gf-5 atomic dataset using the simple and detailed ionisation treatment, respectively. The blue line was 
calculated with the big_gf-4 atomic dataset and the detailed ionisation treatment. The top panels show (from left to right) the normalised 
U,B,V,R,I,J and H passbands of iBessell. »99a) and .Bessell fc Brett. 1,1988,) . 



Thus, with the simple ionisation treatment, the secondary 
maximum is blended in the first peak such that it is not 
clearly evident in the light curves. 

Taking into account all the effects which have been dis- 
cussed in this section we conclude that a detailed treatment 
of ionisation, consistent with the properties of the radiation 
field, has an important influence and is needed to obtain re- 
liable broad-band light curves and spectra after maximum 
light. In particular, we note that the colour evolution ob- 
tained from calculations using a simple LTE description of 
ionisation is subject to strong uncertainties. 



3.2 Influence of atomic data 

As already mentioned, the vast number of line transitions 
associated with the iron group elements are the main con- 
tributors to the opacity. Furthermore the lines provide an 
efficient way to redistribute radiation from the ultraviolet 
and blue pa rt of the spectrum into the red and NIR by 
fluorescence (jPinto fc Eastman! I2OO0I ). This is illustrated in 
Figure |8] which shows, for 20 days after explosion, the wave- 
length at which escaping photon packets were emitted versus 
the wavelength at which the same packets had been last ab- 
sorbed prior to escape (using atomic dataset cd23_gf-5 and 
the detailed ionisation treatment). It shows that many ab- 
sorptions in the UV are followed by emission in the red and 
NIR, but also that the reverse process hap pens too (see dis- 
cussion bv lMazzali|[2000l : ISauer et al.ll2008l ). From the colour 



coding we see that the iron group elements are most effective 
in this redistribution. 

Figure [5] shows the redistribution at 35 days after the 
explosion with escaping photon packets binned into a wave- 
length grid, and indicating how many packets escaped in a 
bin by its grey shade. Thus individual strong lines become 
visible. The dark dots at Aout ~ 8600 A are associated with 
emission in the Ca 11 NIR triplet which results from absorp- 
tion in the Call H and K lines at Ain ^ 4000 A as well as 
from resonance scattering in the NIR triplet itself. In a sim- 
ilar way the Si 11 line at ^ 6355 A shows up as a resonance 
line and in fluorescence with photons absorbed at ~ 4000 A. 
This underlines that a fully detailed treatment of line for- 
mation is needed to make reliable predictions of individual 
spectral features and to a lesser extent also for the broad- 
band light curves. Especially the / band light curve will be 
affected since it has a strong contribution from the Ca 11 NIR 
triplet (see top panels of Figure [5| . 

To simulate this redistribution properly an atomic 
dataset as complete as possible is desirable. Unfortunately 
the simulations become more and more expensive for big- 
ger atomic datasets and the atomic data is not completely 
known. Thus in the following we study the effect of incom- 
plete atomic data by using datasets of different completeness 
to calculate the spectral evolution of W7. 

We started with the line list of iKurucz fc Belli (Il995l ) 
and tried to reproduce the results of the previous section 
using atomic datasets with cuts of -2 and -3 in log (gf) com- 
pared to -5 which was used so far (see Table [1] for details on 
the atomic data). The light curves obtained with the -2 cut 
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Figure 7. UVOIR bolometric and U,B,V,R,I,J,H,K light curves for the W7 explosion model. The green/red lines are the light curves 
obtained with the simple/detailed ionisation treatment and atomic dataset cd23_gf-5. The blue lines are calculated with the detailed 
ionisation treatment and atomic dataset big_gf-4. The small fluctuations e.g. in the U light curve of the simple ionis ation treatment and the 
NIR light curves are due to Monte Carlo noise. For comparison observa tions of SN 2001el |Krisciunas et al.l l20o3) are overplotted as grey 
circles assuming a distance modulus of 31.54 and the reddening law of lCardelli et al" with Av = 0.5 and Ry = 3.1. Furthermore 

W7 light curves obtained with the Monte Carlo radiative transfer code SEDONA (Kasen priv. comm.; UVOIR,U,B,V,R,I,J,H,K dotted 
lines) and the radiation hydrodynamics code STELLA (Sorokina priv. comm.; UVOIR,U,B,V,R,I dashed lines) are shown. 



differ strongly from the previous ones, indicating that lines 
weaker than -2 in log (gf) play an important role. With the 
-3 cut, the light curves in the optical bands are almost iden- 
tical but there are still some minor differences in the NIR 
bands. Going further to a -20 cut does not change the results 
compared to the -5 cut. So the cd23_gf-5 atomic dataset, 
which was used above, is sufflcient to simul ate the redistri- 
bution which is possible within the line list of lKurucz^ fc Belli 
(| 19951 ). 

However, this line list with its meagre total of ~ 5 • 10^ 
lines, contains only a tiny fraction of the millions of lines 
expected of the iron group elements. As many of these lines 
are associated with transitions amongst highly excited lev- 
els, this line list lacks especially lines in the red and NIR 



parts of t he spectrum. It was pointed out by iKasen et al.l 
(2006) and lKasenI (|2006l ) that these missing lines are crucial 
to obtain light curves in good agreement with observations, 
especially in the NIR. To study this, we use atomic dataset 
big_gf-4 (for details see Table [T] and Section I2.3|l which in- 
cludes atomic data for Fell, Felll, Coll, Com, Nill and 
Nilll from more modern theoretical computations, increas- 
ing the number of lines by a factor of ~ 20 compared to the 
cd23_gf-5 atomic dataset. 

The bolometric and broad-band light curves obtained 
for the W7 model using the detailed ionisation treatment 
and atomic dataset big_gf-4 are shown by the blue line in 
Figure [7] UVOIR bolometric and the U, B and V hght 
curves are not strongly affected. In contrast, the enhanced 
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Figure 8. Wavelength redistribution in the last line transition of 
photon packets which escaped at 20 days after the explosion. Each 
point displays the wavelength of an escaping photon packet versus 
the wavelength at which the same packet had been absorbed im- 
mediately before emission. Thereby a single point can represent 
multiple photon packets which have undergone the same absorp- 
tion and emission processes before they escaped. Points on the 
diagonal line are photon packets which last undergone resonance 
scattering. Points above (below) that line have been redistributed 
into the red (blue) by (reverse) fluorescence. The colour coding 
shows by which element an escaping photon was last emitted 
(atomic numbers refer to the colour on top of the number centre) 
thus indicating that the iron group elements are most efficient 
in redistributing flux into the red by fluorescence. The data were 
extracted from our W7 calculation using the detailed ionisation 
treatment and the cd23_gf-5 model atom. 



redistribution of flux from the UV and blue to the red and 
NIR at early times, which can be seen comparing the red 
and blue curves in the top panels of Figure [S] increases the 
first peak in the R, /, J, H and K bands and causes strong 
differences for t < 40 days. As a consequence of this redis- 
tribution, the U and B band light curves using the big_gf-4 
atomic dataset are slightly dimmer between ~ 10 and ~ 40 
days than those using the cd23_gf-5 atomic dataset. At later 
times, the Ught curves are very similar. 

The increase in the magnitude of the primary peak in 
the NIR brings the model light curves to better agreement 
with observations. The distinct secondary maximum which 
is present in the R band light curve using the cd23_gf-5 
atomic dataset changes - as it is observed - to a shoul- 
der using atomic dataset big_gf-4. However it occurs too 
early compared to the observations and especially the H and 
K bands show systematic differences from the observations 
having secondary peaks brighter than their primary peaks. 
This could be due to a further lack of atomic data, which 
might increase the redistribution into the H and K bands. 
However, calculations with the big_gf-3 atomic dataset (see 
Table [TJ, which has fewer lines by a factor of 2, do not 
show significant differences around the first peaks, suggest- 



Figure 9. Wavelength redistribution in the last line transition 
of photon packets which escaped at 35 days after the explosion. 
The escaping photon packets have been binned into a wavelength 
grid. The shade of a grid point indicates how many photon pack- 
ets escaped in this bin (the darker the more packets escaped). 
Points on the diagonal line are photon packets which last under- 
went resonance scattering. Points above (below) that line have 
been redistributed into the red (blue) by (reverse) fluorescence. 
The dark dots at Aout ~ 8600 A and ~ 6300 A show the flux 
redistribution from the Call H and K lines into the Can NIR 
triplet and in Si n respectively. The data were extracted from our 
W7 calculation using the detailed ionisation treatment and the 
cd23_gf-5 model atom. 



ing that very weak lines are not critical in this. It is possible, 
however, that the atomic data set in use is still not complete 
in terms of lines stronger than a -4 cut in \og{gf) - which 
we cannot exclude - or that more iron group elements or 
ions should be inc luded from the comprehensive atomic data 
set (|Kurucj|200^ ). However memory limitations and com- 
putational costs become a severe issue if the atomic data 
are to be significantly expanded and we note that the small 
abundances of the other iron group elements suggest that 
they are less likely to have a major effect. Alternatively, the 
explosion model itself may be at fault. In particular in a 
model where the recombination front from doubly to singly 
ionised material hits the iron-rich inner core later, the sec- 
ondary maxima would be dimmer and in better agreement 
with observations. However, investigating more modern ex- 
plosion models is beyond the scope of this work and will be 
addressed in future studies. 



3.3 Comparison with other codes 

We now compare our bolometric and broad-band light 
curves with those obtained with other radiative transfer 
codes for the W7 model. We start with the light curves 
from the multi-energy group r adiation hydrodynamics code 
STELLA (|Blinnikov et al.l |2006| ) which tre ats the line opac- 
ity in the expansion opacity formalism l|Karp et al ] Il977l : 
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iFriend fc Castodl 19831 : lEastman fc Pintoll 19931 ) . E. Sorokina 
(priv. comm.) provided us with bolometric, U, B, V, R and I 
light curves calculated with an updated version of this code 
using a line hst of ~ 1.6 ■ 10"" lines. In this approach, mat- 
ter is treated in LTE and redistribution is modelled using 
an approximate source function. The light curves are shown 
in Figure [7] aside from having earlier peaks, even in bolo- 
metric light, these light curves are very similar to our light 
curves calculated in the simple ionisation treatment with 
atomic dataset cd23_gf-5. This is not surprising since our 
simple ionisation description is appropriate for LTE and the 
atomic data sets are of comparable size. The obvious differ- 
ence around the first peak in the / band, which is dominated 
by the Ca ll NIR triplet, is most likely due to our more com- 
plete treatment of fluorescence. Compared to our detailed 
ionisation treatment the same remarks as made in Section 
[3TT] apply. 

We also show light curves obtained with the 3D Mont e 
Carlo radiative transfer code SEDONA (|Kasen et al.ll2006l ). 
SEDONA uses the expansion opacity formalism but is capable 
of treating fluorescence with a more sophisticated approx- 
imation for a subset of lines using a downward-branching 
scheme. When a photon is absorbed by a line treated in the 
expansion opacity formalism a two-level atom (TLA) ap- 
proximation is used and the photon undergoes coherent scat- 
tering with probability 1 — e. Otherwise it is absorbed and 
reemitted at a wavelength sampled from the local thermal 
emissivity, thus representing both true absorption and fluo- 
rescence. In principle e is a unique parameter for each line, 
ho wever SEDONA usua lly uses a common value for all lines 
[see lKasen et al" for a discussion! . Matter is treated in 

LTE with the local kinetic temperature being derived from 
balancing the thermal emissivity with the energy deposition 
by 7-ray heating and photon absorption recorded during the 
Monte Carlo simulations. The light curves shown in Figure[7] 
were provided by D. Kasen (priv. comm.) and use a line list 
of ~ 10^ lines treated in the TLA approximation. 

The light curves are in good agreement with our cal- 
culations using the big_gf-4 atomic dataset up to maximum 
light. After maximum light some differences manifest. The 
SEDONA U and B light curves after 50 days decline faster 
than ours - this may be because of recombination to the 
neutral state, which is responsible for the relatively rapid 
decay of the U and B light curves obtained with our simple 
ionisation treatment (see Section|3TT|. Note that, despite us- 
ing LTE conditions to describe the matter state, the SEDONA 
light curves do not fade as quickly after maximum light as 
those obtained with Stella or with our simple ionisation 
treatment (which is equivalent to treating matter in LTE). 
This is likely to be due to the different manner in which the 
kinetic temperature is computed by the different codes. 

The V band light curves are almost identical. In the 
NIR light curves we get remarkably similar results up to ~ 
40 days. Afterwards, the SEDONA light curves stay brighter. 
Since we have seen in Section [3. II that a more complete re- 
combination leads to an increased redistribution of flux into 
the red and NIR by Fe I and Co I this is presumably due to a 
more complete recombination compared to our calculation. 
The strong differences in the / band light curve are most 
likely due to the different treatment of line formation. The 
flux in the / band is dominated by the Call NIR triplet 
for which the SEDONA light curves shown here assume pure 



scattering (e — 0, see the discussion in lKasenll2006l ). while 
we treat the line formation in full detail. 



4 APPLICATION TO AN ELLIPSOIDAL TOY 
MODEL 

In this section we use an ellipsoidal toy model to 
demonstrate the multi-dimensional capabilities of the 
code and to illustrate the basic effects which large- 
scale ejecta asymmetries introduce in spectra and light 
curves. Large-scale asymmetries in the ejecta of SNe 
la are suggested both by observed polarimetry [see 
IWang fc Wheeled (|2008l ) for a review] as well as by 
theoretical explos ion models [e.g. from an off-centre ig- 
nition condition (|Ropke. Wooslev fc H illebrandt' '2 007|), a 
defla gration to detonation transition (Ropkc & Ni emeverl 
l2007f) or in the gravit ationally confined explosion sce- 
nario l|Plewa. Calder fc Lamh,.2004i )] . We stress that the toy 
model we employ here is extremely simple and has a much 
stronger asphericity than the observations or theoretical ar- 
guments suggest. However, it is a useful test case for our 
code which clearly identifies the sense of the effects acting 
on the spectra and light curves. 

4.1 The model 

Taking the total mass and composition of the W7 model and 
assuming homologous expansion, we constructed a simple 
ellipsoidal toy model assuming rotational symmetry about 
the 2:-axis. The maximum velocities along the x- and y- 
axes are set to be only half the maximum velocity («™^^ = 
27 500 km s~^) along the 2:-axis, thus giving an axis-ratio of 
1:1:2. Within the model, ellipsoidal surfaces are taken as 
surfaces of constant density and we adopt a density profile 

p{v) tx exp + (y) /""o^ (19) 

in cylindrical polar coordinates and vo — 2750 km s"'^. We 
assume a stratified composition with three zones. The in- 
nermost zone contains all the iron group material (Sc to Zn) 
and is surrounded by a zone of intermediate mass elements 
(F to Ca). The outermost zone contains the unburned ma- 
terial (IIe,C,N,0). Inside each zone the relative abundances 
are homogeneous and kept fixed at their W7 values. 

We mapped this model to a 64^ grid and followed the 
propagation of 3.2 ■ 10^ energy packets over 100 time steps 
from 2 to 80 days after the explosion, using the cd23_gf-5 
atomic dataset and the detailed ionisation treatment. The 
first 10 time steps have been calculated in the simple ionisa- 
tion treatment to get reliable initial values (see the discus- 
sion in Section r2.2.2|) . To save computational time we make 
use of our initial grey approximation setting TTh.min ~ 15 
and Ngray = 40 (see Section [27lT6| . Spectra and light curves 
were obtained by binning the escaping packets in 10 equal 
solid angle bins centred around the z-axis. 

4.2 Spectral evolution 

In Figure [Tol we compare the 2500 to 10 000 A spectra along 
the major and minor axes with an angle averaged spectrum 
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Figure 10. Spectra for the ellipsoidal toy model at 15, 25 and 
40 days after the explosion (from top to bottom). Blue/red lines 
are for viewing down the major/minor axis. The black line shows 
an angle- averaged spectrum for comparison. 

for 15 (around maximum light in B band), 25 and 40 days 
after explosion. In the absorption trough of the P Cyg fea- 
ture of the Call NIR triplet one clearly sees the different 
velocity extent along the different axes: viewed down the 
major axis, the escaping photons see a velocity field twice 
as large compared to the minor axes such that the absorp- 
tion troughs extend farther into the blue. In general the 
spectra along the minor axes show much sharper features 
than those along the major axis, where the higher velocities 
cause stronger blending. 

Furthermore the total flux observed along the minor 
axis is larger. This can be understood by a simple geo- 
metrical argument. Viewed down the minor axis the cross- 
section surface of our ellipsoid is twice as big as seen down 
the major axis. For an opaque ellipsoid of uniform surface 
brightness, we would expect the same ratio for the flux (i.e. 
AMgeom ~ —0.75 in magnitudes). However, additional ef- 
fects come into play such that the flux difference between 
major and minor axis depends on wavelength and time as 
discussed below. 

4.3 Broad-band light curves 

The broad-band light curves are shown in Figure 1111 As 
expected, the light curves observed along the minor axes 



Table 2. AM = Mminor — A^major fof Selected bands in the 
ellipsoidal model at different times after explosion. 



t 


15d 


22d 


29d 


35d 


42d 


49d 


AM{U) 


-0.93 


-1.01 


-0.92 


-0.80 


-0.70 


-0.87 


AM{V) 


-0.52 


-0.65 


-0.54 


-0.48 


-0.45 


-0.46 


AM{I) 


-0.42 


-0.41 


-0.35 


-0.36 


-0.34 


-0.23 


AM{H) 


-0.38 


-0.24 


-0.22 


-0.27 


-0.27 


-0.19 



Table 3. Peak times and AM15 for selected bands in the ellip- 
soidal model. 



(d) t-LT^ (d) AM^"5^J°' AM\ 



UVOIR 


14.7 


15.8 


-0.6 


-0.9 


u 


11.8 


13.7 


-1.7 


-1.9 


B 


15.2 


15.8 


-1.7 


-1.8 


V 


18.3 


19.0 


-0.9 


-1.1 


i?. 


19.0 


19.0 


-0.5 


-0.6 



are always brighter than those observed along the major 
axis. The viewing-angle effect is always strongest in the bluer 
bands and after maximum light becomes weaker with time 
in V and redder bands (see Table [2]) as the ejecta become 
optically thin in these parts of the spectrum. 

In the U and B bands, which stay optically thick 
throughout our calculation, this effect remains significant 
until the latest times of our simulation and the difference 
AM = Afniinor — Mmajor IS evcu bigger thau one would ex- 
pect from the simple geometrical argument (which would 
imply AMgeom ~ —0.75; see above). This is simply because 
we do not have an opaque ellipsoid of uniform surface bright- 
ness. In fact. Figure [12] - which illustrates the region of last 
emission (RLE) of selected bands at different times - shows 
that the U band RLE is concentrated around the equatorial 
plane which acts to amplify the geometrical effect (as long as 
the ejecta are optically thick, the direction of escape tends to 
be peaked normal to the contours of constant density) . This 
concentration is a consequence of the large number of highly 
optically thick Fe group lines in the UV: U band photons 
are more strongly trapped than photons in other bands and 
therefore tend to preferentially leak out along the equatorial 
plane where the velocity swept out is smallest. 

In contrast, the RLE of V and redder bands do not show 
a strong enhancement around the equatorial plane. The V 
band RLE around maximum light is nearly ellipsoidal, mak- 
ing AAf close to AAfgcom. At later times, the emission be- 
comes more isotropic as the ejecta become optically thinner 
(in Figure [12] the whole ellipsoid contributes to the RLE) 
and AM decreases. For the / and H band, in which AA/ is 
always less than AAfgcom, the ejecta are already becoming 
optically thin around maximum light. 

The H band RLE shows a slight enhancement towards 
the polar regions around maximum light. This is related 
to the reduced U band emission in these regions. The high 
optical depths for blue and UV photons mean that fluores- 
cence redistributes flux into the NIR where optical depths 
are lower such that the radiation can escape. 

In general the NIR light curves are least viewing-angle 
dependent, thus supporting the use of NIR light curves for 
cosmological distance measurements since less intrinsic scat- 
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Figure 11. UVOIR bolometric and U,B,V,R,I,J,H,K light curves for the eUipsoidal toy model. The blue/red lines are the light curves 
obtained along the major/minor axis. The black line shows an angle- averaged light curve for comparison. The filled circles in the //-band 
panel indicate the times for which snapshots of the ionisation state of Fe are shown in Figure 1131 



ter would be expected if geometry effects have any role in 
the observed properties of SNe la. However more detailed 
studies are needed to investigate this. 

The Ught curves observed along the major axis peak 
slightly earlier than those observed along the minor axes (see 
Table Furthermore Ught curves observed along the major 
axis decline more slowly than those observed along the minor 
axes (compare the AMulfl values in Table|3]). A similar effect 
wa s already fo und in a study using a grey version of our code 
bvlsiinl (|2007| ). and we note that it is opposite to the sense 
of the observed light curve width-luminosity relation. 

The plot showing the regions of last emission in Fig- 
ure [12] also traces the stratified composition of our model. 
The U and / bands - which have non-negligible contribu- 



AAfi5 gives the change in magnitude between maximum light 
and 15 days after maximum light. 



tions of Ca due to the Ca ll H and K lines in U and the Ca ll 
NIR triplet in 7 - show significant emission from the zone of 
intermediate mass elements. In contrast, the RLE in the V 
and H bands - which are dominated by Fe group elements 
- is concentrated towards the iron-rich inner core. 



4.4 Secondary maximum in the NIR bands 

As discussed by iKasenI (|2006l ). the secondary maximum in 
the NIR bands forms when the inner iron-rich core recom- 
bines from doubly ionised iron group elements to the singly 
ionised state. With the singly ionised state becoming dom- 
inant, the number of optically thick lines in the blue in- 
creases, thus blocking blue radiation very effectively and 
making the redistribution of fiux from the UV and blue part 
of the spectrum into the red and NIR by fluorescence very 
effective. Here we discuss this effect for the H band light 
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Figure 12. Region of last emission for selected bands ((7, V , I, H from top to bottom) and different times (from left to right). Dark 
regions contribute most to the fiux escaping in the band. The solid lines indicate the three zones of different composition, with Sc-Zn in 
the centre surrounded by F— Ca and finally He,C,N and O in the outer zone. The model is symmetric under rotation about the z-axis. 
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Figure 13. lonisation fractions of Fell, Felll and Feiv (from top to bottom) for selected times (from left to right). The corresponding 
times are marked with filled circles in the H band panel of Figure [TT] Black regions have an ionisation fraction of 1. Lower ionisation 
fractions are shown in lighter reddish colours down to white (0) on a linear scale. 



curve and show the ionisation fractions of Fell-iv (which 
we use as a proxy for all iron group elements) in Figure [13] 
for selected times from near the first peak (~ 10 days after 
explosion) to well after the second peak in the H band (~ 
32 days after explosion, compare Figure [11} . 

Around the first peak the iron-rich core consists pre- 
dominantly of Feiv. In the very centre there is still a con- 
tribution of Fe V (not shown in Figure I13p and at the outer 
edge there is a small region which has already recombined 
to Feiii. As time passes the ejecta cool and the Feiv region 
recedes. At the same time the Felll region extends further 
in and a contribution of Fe ll starts to build up around the 
outsides. Recombination occurs first at the poles since they 
see the least amount of ionising radiation (compare Figure 
I12p . The secondary maximum in the H band occurs at ~ 32 
days when the Fen region first forms a complete elliptical 



ring (Figure [T3)l . This is the point at which the opportunity 
for fiuorescence is maximal. Note that although the ionisa- 
tion fractions vary with position, the NIR bands are suf- 
ficiently optically thin that the NIR fiux is approximately 
angle-independent at these epochs (this is apparent from 
Figure [121 where the full ellipsoid is bright). At late times 
the recombination front moves inwards as the light curves 
fade until only a small Fe ill core remains at the end of our 
simulation. 

Comparing the two different lines-of-sight in Figure [TT] 
we see that the local minimum between the first and sec- 
ondary maximum in the H band light curve along the po- 
lar axis occurs slightly earlier (~ 15 days) than along the 
equatorial axes (~ 18 days). This follows from the earlier 
recombination from Fe ill to Fe ll in polar lobes which gives 
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an earlier increase in the redistribution of flux from the UV 
and blue into the red and NIR. 



5 CONCLUSIONS 



Based on the approach of ILucvI (|2002| . |2003| ) we extended 
the sre y time-dep endent 3D Monte Carlo radiative transfer 
code of fsiinl (|2007l ) to a non-grey opacity treatment. The new 
code, ARTIS, treats the 7-deposition and spectrum forma- 
tion in detail and solves the ionisation balance together with 
the thermal balance equation consistently with the radiation 
field. Line formation is treated in the Sobo l ev approxim a- 
tion using the macro- atom approach of lLucvl (|2002l . l2003h to 
model atomic physics in detail. This allows parameter-free 
radiative transfer simulations for 3D hydro models with a 
maximum of predictive power to be made. 

We applied this code to th e well-studied one- 
dimensional deflagration model W7 l|Nomoto et al.l Il984l ) 
as a test case and found good agreem ent with both ear - 
lier wor k using different codes [sedona dKasen et al. 



STELLA iBlinnikov fc Sorokina| 2002|: Blinniko v et al 



200e) 



2006D ] 



and observations fSN 1 994D ijPatatet al. 1996i); SN 2001el 
ijKrisciunas et al.ll2003l )]. Concerning the NIR light curves, 
we confirm the importance of line fluorescence in modelling 
these bands and, particularly their strong dependence on 
the size of the atomic data set in use. Rather than adopt- 
ing LTE, we assume photoionisation equilibrium and solve 
the ionisation balance equations using rates consistent with 
the radiation field. We showed that this leads to significant 
differences in the ionisation structure after maximum light, 
which strongly affect spectra and light curves. The macro- 
atom formalism allows us to avoid introducing a parame- 
terised treatment of line fluorescence, such that we account 
for the differing contributions of resonance scattering and 
fluorescence in different lines (see Figure |9}. 

To demonstrate the multi-dimensional capabilities of 
the code we have presented calculations for an ellipsoidal 
toy mode l as an exa mple. As expected [e.g. from the grey 
study by [si^ (|2007l ) for a similar model], light curves ob- 
served along the minor axes are brighter than those observed 
along the major axis. The sensitivity decreases with time, 
as the ejecta become less optically thick, and from blue to 
red wavelengths. If line-of-sight effects due to asymmetric 
explosions have a significant influen ce on the scatter around 
the Phillips relation of SNe la (Phillipq 199:^) ■ this could 
explain why the NIR light curves seem to be more homo- 
geneous and thus most promising for cosmological distance 
measurements. However more detailed studies are needed to 
investigate this fully for realistic models. 

In future work we will extend our method to include 
a treatment of non-thermal excitation and ionisation which 
we currently do not treat in detail but which can be impor- 
tant (for a discussion see lBaron et al.lll996h . This will also 
allow us to extend the simulations to later epochs (nebu- 
lar phase) which are strongly affected by these processes. 
Eventually time-dependent terms in the statistical equilib- 
rium equations could be taken into account and polarisation 
could be treated, thereby extending the range of data with 
which we can make meaningful comparisons. 
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